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High-energy theory for close Randall Sundrum branes 
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We obtain an effective theory for the radion dynamics of the two-brane Randall Sundrum model, 
correct to all orders in brane velocity in the limit of close separation, which is of interest for studying 
brane collisions and early Universe cosmology. Obtained via a recursive solution of the Bulk equation 
of motions, the resulting theory represents a simple extension of the corresponding low-energy 
effective theory to the high energy regime. The four-dimensional low-energy theory is indeed not 
valid when corrections at second order in velocity are considered. This extension has the remarkable 
property of including only second derivatives and powers of first order derivatives. This important 
feature makes the theory particularly easy to solve. We then extend the theory by introducing a 
potential and detuning the branes. 
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I. INTRODUCTION 



Motivated by advances in string theory (in particular 
heterotic M-theory), there has recently been consider- 
able interest in models where spacetime is effectively five- 
dimensional. In these models, matter fields are confined 
to three-branes, membrane-like surfaces embedded in the 
higher dimensional space, while gravity (and other bulk 
fields) can propagate in the whole of spacetime Q . The 
cosmological consequences of these scenarios have been 
widely studied (for recent reviews see @,O1>0)- The sim- 
plest such model is that of Randall and Sundrum 0, [|J 
where the bulk is assumed to be empty apart from a cos- 
mological constant. This is a toy model through which 
braneworld ideas can be tested and, despite its simplicity, 
is rich in new physics. 

In the low-energy limit, an effective four-dimensional 
theory can be derived on the branes BUSES This 
theory predicts that near the collision, the Hubble con- 
stant on each brane is related to the proper contraction or 
expansion velocity of the fifth dimension d by H± = jp- , 

radius related to 



where L is the Anti-de Sitter (AdS) 
the bulk cosmological constant. However, an exact cal- 
culation gives the result H± = tanh 2 (ij- As ex- 
pected, this agrees with the four-dimensional effective 
theory only at low velocities. To lowest order in veloc- 
ity, the low-energy limit gives an accurate result for the 
brane geometries. The aim of this work is to go beyond 
the low-energy limit and to develop a covariant formal- 
ism which describes exactly these velocity corrections in 
the small distance limit d<L. 
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Braneworld cosmology offer a entire new set of pos- 
sibilities for the production of the large scale structure 
and in particular it allows scenarios with high energies 
or for which the branes are moving apart at velocities 
which could be large HH [13, 113 E3, EH [IE E3, EE 
I20L EH I22I . We may therefore address the issue of 
the importance of the high-order terms in these regimes. 
In scenarios such as steep inflation [22I I23I ] , the energy 
scales are important and facilitate the end of inflation. 
In the Cyclic Universe 0, US as another example, it is 
suggested that terms of fourth order in velocities should 
be considered in order to obtain a scale invariant spec- 
trum after the bounce j^S- To work with such models, 
it seems important to understand the behaviour of these 
high-energy corrections and in particular to understand 
their consequences for the production of the large scale 
structure, as much in the context of an expanding as in 
a contracting Universe. 

However, in order find a solvable theory when the low- 
energy constraint is relaxed, we need to work in another 
limit. The special limit in which we choose to work in- 
stead is the close brane limit where the distance d be- 
tween the branes is small compared with the AdS radius 
L: d <C L. This limit is interesting in cosmology for two 
main reasons. The first one is that there has been con- 
siderable interest in the interpretation of the Big Bang 
as a brane collision IH IM 13, El El El US ■ There- 
fore, at the beginning of our universe, there will be a 
regime where the distance between the branes is very 
small. Such a limit seems therefore to be relevant to 
study the production of the large scale structure. 
The second reason is that it is in the limit where the 
distance between the branes is small that we expect 
the geometry on the brane to be well-described by a 
four-dimensional theory. Indeed, when the branes are 
far apart, the bulk degrees of freedom (leading to the 
Kaluza-Klein corrections on the branes) are more eas- 
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ily excited, making the theory on the branes non-local 

Following these arguments, in the small distance limit, 
we therefore expect the branes to be well described by a 
four-dimensional theory which will be valid for any en- 
ergy scale and will include higher-order derivatives. 
To work in this re gim e, we follow part of the idea of 
Shiromizu et.al. in [34] although our final result will be 
different. The main idea of this paper is to express the 
extrinsic curvature on the negative-tension brane as a 
Taylor expansion in terms of the extrinsic curvature on 
the positive-tension brane and its derivatives along the 
normal direction. We can then use the five-dimensional 
Einstein equations to express any second (or higher) or- 
der normal derivatives of the extrinsic curvature in terms 
of first normal derivatives and derivatives along the four 
transverse directions. Since the normal derivative of 
the extrinsic curvature on the positive-tension brane is 
known up to the induced Weyl tensor on that brane, this 
gives a formal equation for the induced Weyl tensor on 
the positive-tension brane, which is the only unknown 
information on the brane. 

Although formally correct, this process can be compli- 
cated, in particular if matter is introduced on the branes. 
In order to understand this expansion we will therefore 
focus on the case of empty branes for which the extrinsic 
curvature on the branes is proportional to the metric. As 
it will appear later, the formalism becomes very simple in 
that context. As already mentioned, the second condition 
we impose is the close brane regime. More concretely, by 
imposing this condition we will keep only the terms at 
leading orders in the distance d between the branes at 
each order in the Taylor expansion. To do this, we will 
work in a recursive way. This will enable us to work out 
the exact form for the Weyl tensor on the positive-tension 
brane to leading order in d. 

As we might expect the exact expression for the Weyl 
tensor in the close brane limit has some higher derivative 
correction terms which are not present in the low-energy 
four-dimensional effective theory. However we may check 
that, at low velocities, this expression for the Weyl tensor 
is consistent with the one derived from the low-energy ef- 
fective theory in the small d limit. Furthermore, as a non 
trivial check we may verify that our result gives precisely 
the right result to leading order in d for the background 
solution where it is possible to solve the five-dimensional 
geometry exactly. However the prescription is completely 
covariant and will allow us to study perturbations with- 
out needing to solve the full five-dimensional equations. 

Having an exact expression for the induced Weyl ten- 
sor in the close-brane limit provides us with a unique 
modified Einstein equation on the brane describing the 
coupling of a scalar field (the radion) to the brane ge- 
ometry. The equation of motion of the scalar field is 
then given by the requirement that the Weyl tensor be 
traceless. One of the most promising routes is the possi- 
bility of interpreting this scalar field as playing the role 
of the inflaton scalar field on the positive-tension brane 



[H IH 0i |H HJ. During inflation both branes 
could be moving apart. When the negative-tension brane 
moves towards infinity, its effect on the positive-tension 
brane becomes negligible which means that the scalar 
field decouples, giving an explanation of why such a field 
is not seen on the brane at the present time. In order 
to interpret the radion as a candidate for the inflaton a 
potential for the radion has to be introduced. So far the 
exact origin of that potential is not fully understood, but 
for the purpose of this study we shall introduce a po- 
tential by hand assuming some five-dimensional effects 
might explain its presence. We will therefore study how 
the theory we derived on the brane can be modified in 
order to include such a potential and how the equation 
of motion for the scalar field has to be modified in a con- 
sistent way in order to respect the conservation of energy 
and to be consistent with the low-energy theory. In order 
to do so we will review how the introduction of the po- 
tential is usually performed in the low-energy limit and 
then proceed in a similar way for the close-brane theory. 

The paper is organised as follows. In Section [U] we 
review the Randall Sundrum model and derive the ex- 
act Friedmann equations on the branes for the back- 
ground. We then give an overview of the low-energy four- 
dimensional effective theory and compare its predictions 
in the background with the exact solution. Since they 
disagree beyond the leading term in velocity, in Section 
IIIII we derive a covariant expression for the Weyl tensor 
on each brane in the close brane limit. First we start 
with a toy model in order to understand the procedure. 
Then we present the five-dimensional model and derive 
the Weyl tensor on each brane using the small d approx- 
imation d <C L. We use this result in Section IIVI in 
order to find the exact theory on the branes, valid in the 
close brane limit. We check that this theory is consistent 
with the conservation of energy. Finally in Section Ivlwc 
consider extension of this theory, first by introducing a 
potential and then by detuning the brane tensions. We 
will then discuss the implications of our results in Section 

eh 

In Appendix [S] we present the technical details that 
allows us to derive the normal derivative of the extrinsic 
curvature in order to obtain the Weyl tensor on each 
brane. First we proceed in a special gauge where we 
suppose that the g vy component of the metric can be 
taken to be independent of the fifth dimension. Then, we 
show that the same procedure is valid when this condition 
is relaxed and we obtain the same result. 

In Appendix[Bl we present the details of the derivation 
of the divergence of the stress-energy tensor and show 
how the equation of motion for the scalar field has to be 
modified when a potential is introduced. 

Finally, for completeness, we derive the evolution equa- 
tion for the Weyl tensor in Appendix IU1 
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II. BACKGROUND BEHAVIOUR 



so the brane velocities must satisfy: 



In what follows, we shall be interested in the two brane 
Randall Sundrum model as a specific simple example of 
braneworld cosmologies. In that model, the spacetime is 
five dimensional, with a compact extra dimension having 
the topology of an S1/Z2 orbifold. The stress energy of 
the bulk is assumed to be from a pure negative cosmo- 

logical constant |A| = -jj, k — = -hj, with M n the 

5 4 

n-dimensional Planck mass. For simplicity, we set k — 1 
for the rest of this work. There are two boundary branes 
located at the fixed points of the Z2 symmetry on which 
gauge and matter fields are confined. Apart from where 
stated otherwise, we assume, for simplicity, that the ten- 
sion on each brane are fine-tuned to their canonical value 
X± = ±4 and we don't include any other kind of matter 
on the branes. 

In this paper we use the index conventions that Greek 
indices are four dimensional, labeling the transverse x* 1 
directions, while Roman indices are fully five dimen- 
sional. We will denote by a "dot" any derivative with 
respect of the proper time unless specified otherwise. 



A. Five-dimensional Theory 

1. Frame where the bulk is static 

In the fine-tuned case with cosmological symmetry (i.e. 
the three spatial directions are assumed to be homoge- 
neous and isotropic) one can make considerable exact 
analytical progress. We work in the frame where the 
bulk is static, which exists as a consequence of Birkhoff 's 
theorem. In that frame, the most general geometry is 
Schwarzschild-Anti-de-Sitter (SAdS) with the parameter 
C associated with the Black Hole mass ^1] : 

ds 2 = dr 2 -n 2 (y)dT 2 + a 2 (r)dx 2 (1) 

with a 2 {Y) = e- 2Y ' L + - e 2Y ' L 

4 

n 2 (Y) = L 2 a'{Y) 2 =a 2 -^, 

where we are assuming flat spatial sections for simplic- 
ity. In this frame, the branes are not static but have loci 
Y = Y±(T). The Israel matching conditions jlj, associ- 
ated with the Z2-symmetry, impose the condition on the 
extrinsic curvature: 



K^{Y = Y ± ) 



1 



g^(Y = Y±). 



(2) 



On the brane, the spatial components of the extrinsic 
curvature are: 



Ku<y±) = 1 



n\ ay) 



a(Y) 



9ij 



(3) 



n 



1 



(4) 



where a±(T) = a(Y — Y±(T)) is the value of the scale 
factor on each brane, and similarly for n±(T). The in- 
duced line clement on the branes can be read off as 



ds| = -{n 2 ± -Yl)dT 2 +4dx 2 



-Ai 



dx 2 ), 



defining the conformal time on each brane: 



At± 



n 



AT 



'-AT. 



(5) 



(6) 



In terms of the conformal time, the scale factor on each 
brane evolves with constant velocity: 

2 



/do±\ 2 = - 2 /dT\ 2 



da(Y) 



± 



Vdr±y AY 



y=y± 



leading to the Hubble parameter on each brane: 



H 2 ± 



1 

4 
1 

J? 



L 2 a 4 ± 



(8) 



This result is a direct consequence of the projected Ein- 
stein equations on the brane (see Section Till C|l . which, 
for empty Friedmann Robertson Walkers (FRW) branes, 
gives the induced Ricci scalars as 



R{±) = _6_d2a± 
dr^ 



0. 



(9) 



We may point out that in JHJ, the Hubble parameter on 
each brane is bounded by L 2 H± < 1, with equality for 
the negative-tension brane when this one touches the five- 
dimensional Black Hole horizon. This study is beyond 
the purpose of this paper and we consider both branes to 
be far away from the five dimensional singularity. 

We now consider the specific case where the branes are 
moving apart after a collision at t = and we denote by 
do the value of the scale factor when the branes coincide 
(the situation where the branes are moving towards each 
other before the collision is completely analogous). 

In contrast to the exact treatment above we now only 
consider the system in the limit of small brane separation, 
i.e. a+ w a_ 1=3 ao, n+ 1=3 n_ 1=3 no and Y+ « — Y-. Using 
(0J), to linear order we have: 

Y+(T) = Yo-v(T-T )(l + o(^^X) (10) 



L 

Y_(T) = y + v(T - T ) ( 1 + O ( ) ) (11) 



L 



Y=Y± 




(12) 
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where the collision happens at T = To and Y = Yq. 
When the branes are close, the proper distance between 
the branes goes as d ~ (T — Tq), so the corrections are of 
order O (2-^) = O (f). This is the limit in which we 
will work all through this paper. 



2. Frame where the branes are quasi-static 

So far we have worked in the frame where the bulk was 
static. However if we are interested in the proper distance 
between the branes, it is more intuitive to derive it from 
the frame where the branes are static. Such a frame is 
in general complicated, however, for the purpose of this 
study, it is enough to consider the frame in which the 
branes are "quasi-static" , or static to leading order in d/L 
or in (T — Tq)/L. In order to work in such a frame, we 
may perform the gauge transformation (Y, T) — > (y,t): 



T = Tq + — cosh f (2y - 1) tanh 1 

no V V"o 



Y = Y + t sinh (2y - 1) tanh 



(13) 
(14) 



We can indeed check that at y = 0, Y = Yq — v(T — To) 
and at y = 1, Y = Yq + v(T — To). In this new frame, 
the branes are static to leading order in d and located at 
y — and y — 1. In terms of the new coordinates y and 
t the bulk geometry in the limit of close separation is: 

ds 2 = A(tfdy 2 -dt 2 + a 2 dx 2 , (15) 



with Ait) = 2nanhT 1 ( — 
\no 



(16) 



In this limit, the induced metric on the branes is then 
ds 2 - = — di 2 + a(t) 2 dx 2 , where t is the proper time. When 
d <C L, the proper distance between the branes is given 
by: 



d = 



A{t) dy 



2 t tanh 



no 



(17) 



so that, to leading order, the expansion of the fifth di- 
mension with respect to the proper time t is 



d = 2tanh~ 1 ( — ), 
no 



(18) 



where v is given in (|12fl . At the collision, v is related to 
the Hubble parameter (JHJ by: 



2 2 



l-^ = L 2 H 2 ± (r = 0), 



(19) 



This result is derived directly from the full five- 
dimensional equations subject only to the assumption 
that the branes are close. In JHJl , we have already pointed 
out that the Hubble parameter was bounded L 2 H\ < 1 
which is consistent with (|20|l . 

We can compare this result with the analogous rela- 
tion derived from the standard low-energy effective the- 
ory which we shall briefly review in the next Section. In 
Section Hvl we will then compare this result with the one 
obtained from the small-c? theory which we derived in this 
paper and we shall see that they agree perfectly. 



B. Four-dimensional effective theory at low energy 

1. Low-energy effective theory 

At low energies the system is well-described by a four- 
dimensional effective theory 0, 0, I n this limit, 
braneworlds behave like conventional scalar-tensor the- 
ory of gravity where the bulk effects are represented by a 
single scalar field. In the low-energy approximation the 
induced metrics on each brane are conformally related: 



(21) 



and the field ^ is related to the proper distance d between 
the branes by ^ = e~ d / L . The system is closed and 
described by the modified Einstein equation: 



(22) 



= ^ 2 G { +) + 4d^d v y - 2^D^D V ^ (23) 



2^-{d^) 2 )g£K 



where all covariant derivatives are taken with respect to 
gffi- This is in complete agreement with the Gauss- 
Codacci equations in the low-energy limit. The equation 
of motion for the scalar field is given by the requirement 
that Efu/ is traceless: 



= 0. 



(24) 



In terms of the proper distance between the branes, the 
low-energy effective theory becomes: 



2 f 2 



-i (d,dd u d- l -{dd) 2 gW 



Dd 



(ddy 



(25) 



(26) 



giving the relation between the Hubble parameter and 
the d: 



2. Background behaviour of the effective theory 



±— tanh 
L 



for d <C L (20) ^ e mav now examine the behaviour of this effective 

theory in the special case of a flat FRW background in 
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order to compare the results with those of Section Til Al 
We consider the special case where the positive-tension 
brane expands. Since the Ricci scalar vanishes, the scale 
factor still satisfies: 



d 2 a. 



= 0, 



dr 2 



(27) 
(28) 



where ao is the value of the scale factor at the collision 
and i>+ is an arbitrary constant. Comparison with the 
five-dimensional result (JJJ would identify v + as yC/L 
but, considering only the four-dimensional effective the- 
ory, its value is not determined. The Friedmann equation 
obtained from (|23[1 is: 

(1 - * 2 ) Hi - 2f*iJ + - * 2 = 0. (29) 

This is a quadratic equation with the two solutions: 



d 



(30) 



The two possible signs correspond to the branes moving 
either in the same (— ) or opposite (+) directions, which 
can be seen as follows. From the conformal relation i|21|) 
(i.e. a_ = ^a+), the Hubble parameter on the negative- 
tension brane can be written in terms of H + as: 



a_ d 
a_ L 
TH+{l + 0{d/L)) 



+ H+ = Te d/L H+ 



(31) 



So a — sign in '>()[! corresponds to the situation where 
~ H + for which the branes are moving in the same 
direction. In that case, near the collision, 



d w —H + d %-d 



(32) 



so the branes take an infinite time to collide. Instead, we 
are interested in the situation where the collision happens 
at finite time and the branes are moving in opposite di- 
rection H- ~ —H +) corresponding to a + sign in (|30|l . 
In this case, 



H ± =± 



for d < L, 



which implies 



2Lv+ 



(33) 



(34) 



In that case we notice that d is of order (d/L) . We 
may compare this result with 12Ufl , which is correct to all 
orders in d for small d/L. As expected, the low-energy 
theory reproduces only the leading term. The effective 
theory is therefore only valid to leading order in veloc- 
ities (even for the background), but for any value of d. 



We may point out that in this result both LH± and d 
appear to be unbounded. This is only true in this low- 
energy regime for which the restrictions are very strong 
LH± ~ d -C 1. As seen in (|20|) when LH± ~ 1 some new 
restrictions have to be imposed. 

Since the low-energy effective theory only predicts the 
leading order of velocity, it might be possible to go be- 
yond the low energy restriction and find a theory which 
would be valid to all order in velocities. In order to de- 
rive such a theory, we will work in a regime where the 
branes are close to each other. In the rest of this paper 
we show the existence and consistency of such a theory 
which successfully reproduces (|20|l and agrees with the 
low-energy theory in the regime of small separation and 
velocity in which they are both valid. 



III. CO VARIANT APPROACH IN SMALL d 
LIMIT 



A. Toy model 

In order to understand the procedure we will follow to 
work out the theory on the brane in the close brane limit, 
we first examine the following one-dimensional example. 
We consider a second order differential equation: 

f"(y) = U(y)f'(y) + V(y)f(y) + W(y), (35) 

where U, V and W are some known function of y. 
We assume that the function / is known at y = and at 
y = 1 and we wish to find the value /'(0) of its deriva- 
tive at y = 0. One way to do it would be to solve the 
differential equation exactly with the two boundary con- 
ditions for /. Once f(y) is known for any y, we can infer 
/'(0). But by doing so we extracted more information 
that we actually wanted; this method would be equiv- 
alent to solving the five-dimensional Einstein equations 
exactly in order to obtain the induced geometry on the 
brane. Although this is in theory possible it would be 
very hard to do. Instead we will summarise in this ex- 
ample the method used by [34J. The idea is not to solve 
the differential equation exactly but to differentiate it in 
order to use it in the Taylor expansion: 



/(y = i) = EV w 



(0). 



(36) 



n>0 



By differentiating equation (|35|l . we can find an expres- 
sion for f( n '(y): 

= U n ^(y)f(y) + V n - 1 (y)f(y) + W B _i(y),(37) 
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where U n , V n and W n may be found in a recursive way: 

U n+1 (y) = K(y) + V n (y)+U(y)U n (y), 
U 1 = U(y),U = l,U- 1 = 

v n+ i(y) = v£(w) + *W«(y), 

V 1 = V{y),V = 0,V- 1 = l 
W n+1 (y) = <(y) + W%)£/„(y), 

W x = TU = 0, W-! = 0. 

Using these expressions, we may write /'(0) as: 

1 



/'(0) 



/(I) 



(38) 



Knowing £/„, V^j and W n in a recursive way, we can per- 
form the sums and find an exact expression for /'(0). 
This will be very similar to the method we will use to 
find the induced Weyl tensor on the brane. Although 
the extrinsic curvature is known on the branes, its normal 
derivative (which involves the Weyl tensor) is not. We 
can derive, however, a second order differential equation 
for the extrinsic curvature which allows us to calculate 
this derivative in the same way as l|38|) . 

However, already in this linear one dimensional prob- 
lem, the recursive relations are non-trivial. The five- 
dimensional problem is even harder since, unsurprisingly, 
the equations are non linear and formidably complicated. 
However, if we keep only the leading terms in d/L, we 
find that the second order differential equation for is 
linear and, remarkably, the Taylor series corresponding 
to (|35jl becomes tractable. We therefore end up with an 
expression, correct to leading order in d/L, for the nor- 
mal derivative of the extrinsic curvature on the positive- 
tension brane which enables us to write down the Einstein 
equations and obtain the close brane limit of the exact 
theory on the brane. 



B. Regime of validity 

From now on we will work in a regime where the 
branes are very close d <C L. As already mentioned, 
there are two type of solutions for the background, de- 
pending whether the branes are moving in the same or 
in opposite directions. If the branes are moving in the 
same direction, we have seen in (|32[1 that, for the back- 
ground, d ~ d/L. It is subject to this assumption that 
|34j was derived. In this regime, we may check that to 
leading order in d, they recover the low-energy effective 
theory. Therefore, for this regime, the low-energy effec- 
tive theory is valid to all order in velocities (at small d). 
However from l|34|l we see that this is not valid when the 
branes are moving in opposite directions. In that case, 
d ~ ( j\ r-j 1. This will be the solution we will be in- 
terested in for the rest of this work and we will assume 



the relation d^d ~ (jA . Although this is strictly true 
only for the background, if we work with perturbations 
in a comoving gauge, this relation will still hold. In that 
gauge, d^d ~ (4) will still be true covariantly. 

Furthermore, for adiabatic perturbations, the pertur- 
bations behave the same way as the most general back- 
ground solution. For the background there are two kinds 
of solutions, one where the branes move in the same di- 
rection for which d^d ~ i and one where the collision 

happens at finite time, for which d^d ~ This is 

true for adiabatic perturbations as well. The adiabatic 
perturbations will follow a similar evolution to one of 
the two background solutions or a superposition of them. 
Since the low-energy effective theory reproduces correctly 
one type of solution we may focus on perturbations that 
follow the other kind of behaviour for which d^d ~ f-jM . 

At the level of perturbations, one might think that 
this procedure might break down since the perturbations 
diverge. But this divergence is actually logarithmic in d 
and is therefore negligible compared to the terms in 1/d 
that we will find in the theory (|69|l . Compared to L/d, 
\og{d/L) ~ (d/L)°. In 3jj, it is actually shown that in 
the right gauge, the perturbations remain "small" going 
towards the bounce. 



C. Gauss Codacci formalism in the frame where 
the branes are static 

1. Formalism 

In this subsection we will follow the formalism of [34| 
We choose coordinates where the metric is of the form 



ds 2 



g a bdx a dx b 
A 2 {y,x)dy 2 



q^(y,x)dx fl dx 1 ' , (39) 



with x 11 the transverse coordinates and y parameterising 
the extra dimension. In this frame, the branes located at 
the fixed positions y — and y = 1. q^ v {y,x) is the in- 
duced metric of a y — y — const hypersurface. In partic- 
ular, giV(x) = q^iy = 0,x) and 5^ (at) = q^{y = l,x) 
are the induced metrics on both branes. 
The Einstein tensor on a y = const hypersurface will be 
written simply as G^ v , whereas the fully five-dimensional 
tensor will be denoted by ^G^ v . Four-dimensional 
quantities may be expressed in terms of five-dimensional 
quantities by means of the Gauss-Codacci formalism 
[44. |45L Hg • Using the bulk Einstein equations 

{5) G ab = —^g ab , (40) 

the modified four-dimensional Einstein equation is: 

G>„ = A g/ „, + KK^ - K^KS (41) 



1 
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where E^ v is the Electric part of the five-dimensional 
Weyl tensor (see Appendix IO. All indices are raised and 
lowered with respect to the four-dimensional metric q^ v . 
For a Z2-symmetric brane, the extrinsic curvature K^ v 
can be uniquely determined on the brane using the Israel 
Matching conditions, which reduce to: 

KM = 0)=K£( V = 1) = -- St. (42) 

Substituting this into (|41|l gives the projected Einstein 
equation on the branes: 



^ (IV 



_£(±) 



(43) 



where G$ = G^ u (y = 0) and G^' = G^ u {y = 1) are 
the induced Einstein tensors on the two branes. The aim 
of this work is to find E^' exactly in the close brane 
approximation. As shown in 0, [45| the Weyl tensor 
can be expressed in terms of the normal derivative of the 
extrinsic curvature: 



<-) _ 



1 D^D A 1 

- d v K - K»K« + —6», (44) 



where represents covariant derivative with respect to 
Q)j,v We may define the proper distance d between the 
positive-tension brane and the hypersurface located at 
y = const by d(y,x) = L A(y' ,x)dy' . The proper dis- 
tance between the branes is defined as d(x) = d(y = l,x). 

The Lie derivative L„ = = -kd y is the derivative 
along the normal vector of any y — const hypersurface. 
In particular the extrinsic curvature is the derivative 
along the normal vector of the induced metric on such 
a hypersurface: K^ v (y,x) = ^d^q lill (y, x). 

In what follows we will use a similar procedure to the 
toy model. Knowing the extrinsic curvature on both 
branes, we write an expression for its derivative using 
the Taylor expansion: 



K{y = 1) 



E 



y=0 



(45) 



In what follows, we will use the notation: Q' = d y Q and 
Q(ri) = dyQ for any quantity Q(y,x) carrying indices 
only along the directions of the y = const hypersurface. 

To start with we will consider A to be independent of 
y. This is indeed the case for the background geometry 
of the close-brane limit (|15|l . We will see in Appendix IA 21 
that this assumption does not affect the final answer. In 
that case d{y,x) = yA{x), and in particular the proper 
distance between the branes is: d(x) = A(x). We can 
now rewrite the expression (|44|l in the form: 

= -dE£ - D^Dyd - dK£K™ + ±8*. (46) 

In order to find the expression for K t ^ in (|45l) we need 
the derivative of the Weyl tensor and of the Christoffel 
Symbol. The last one is given by: 



pa/ 

1 fJ,V 



D^dK°) + D„(dK°;)-D a (dK^). (47) 



For the Weyl tensor derivative, we may use the result 
derived in Appendix O 

£7f =d(2K2E£ - \kE» - \k^E^ (48) 

-L D * [d 2 D»K av + d 2 D v K£ - 2d 2 D a K»\ . 

where (K 3 ) 1 are some cubic terms in the traceless part 
of the extrinsic curvature which exact form will not be 
relevant for the purpose of this study (they vanish at 
y = and y = 1). Both these relations l|47|) and l|48|l are 
valid for any y. On the brane they simplify considerably: 



r#(o) 

£?t'(0) 



^(d^ + d^-d^ 



Ad 



(49) 
(50) 



The important point is that the system is now closed. 
Writing E^„ back in terms of the extrinsic curvature and 
its derivative, we obtain a second order differential equa- 
tion for the extrinsic curvature which is non-linear but 
which involves four-dimensional quantities only. We may 
therefore apply the procedure of Section fill Al In order 
to do so we will assume the branes to be close and keep 
only the leading order in d in the expansion. 



2. Small d approximation 

In what follow we denote by °K^ n \x) the leading or- 
der in d/L of Kt,^ n \y = 0, x), symbolically, 



K (n) = K (n) ^ + 0(d/L)) . 



(51) 



Using the small distance approximation, we rewrite the 
expansion l|45|l as: 



£ 1 0^(n) =0 



(52) 



where the n = term cancel since the extrinsic curvature 
is the same on both brane l|42|) . All the terms will be kept 
in the sum only if they are all of same order. We shall 
see in what follows that this is indeed the case. In order 
to do so, we will work in a recursive way. 



3. Example for the n=l and n~2 case 

We first concentrate on the n — 1 and n = 2 cases in 
detail in order to gain insight; the technicalities of the 
general n case are left for Appendix lAl 
For n = 1, 



KZ'( y = 0)=-dE2-D»D v d\ v=0 . 



(53) 
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Since d a dpd ~ d a d ~ d° the second term goes as d°. In 
the effective theory, on the brane, E£ ~ (using l|2SJ). 
Although we have argued that at high energy the effective 
theory does not give the exact expression for the Weyl 
tensor, we have seen that (at least for the background) 
the behaviour is the same, differing only in corrections at 
higher order in the velocity. In particular should go as 
d^ 1 at high energies as well (we will check later that this 
is indeed the case). We therefore have K%'(y = 0) ~ d°. 
For the second derivative we have: 



Kt"{y) = -dEt' -q^'DpD v d 
+q" (, Tftd a d-dd v (K£KS). 



(54) 



On the positive-tension brane we may compare these 
terms with the ones from K^'(y = 0): 



terms from K" 
dEf = fE? 

2d ; 



terms from K' 



q»P'DpD v d=2gDi*D u d < D»D u d, 



(55) 



dd v (K£K?) 



2d 



where we used q^' = dd s q^ = -2dK^. When d < L, 
the only term which is not negligible in comparison to 
the terms present in K'%(y = 0) is the one coming from 
the derivative of the Christoffel symbol: 

q^T^(0)d a d = -j(2d»dd„d-(dd) 2 S^ (56) 

This last term will give a contribution of the same order 
as K$'{y = Q). 



4- General n case 

In the previous specific case, we had ~ d° and 

°K"( 2 ) d°. In Appendix lAl we will show that the same is 
true for any n. In particular we will show that the leading 



contribution in Kt '■ (y = 0) comes from the (ri — 1) 
derivative of this Christoffel symbol: 

OK*™ = <f0T$ n - l \o) d a d, (57) 

and the leading contribution from the (n— l) th derivative 
of the Christoffel symbol comes from the term: 

r^-^CO) = d, °K a v {n - 2) +d, u Kf n ~ 2) (58) 
a j\p v 

Using these two results, we therefore have: 

°K* W = d'" Q K a v {n ~ 2) + d, v °K a ^ (n ~ 2) (59) 
-d> a °K^ n - 2 l 

In the Taylor expansion, the term Kt therefore 
comes in with the same contribution as the term 
°Ku^ n ~ 2 \ so all terms have to be considered. 



fch 



5. Expression for the Weyl tensor 

We now use these results in the Taylor expansion. 
First, we may define the operator O such that: 



OZ£ = [d'^ZI + d^Z a » - & a z%\ d, a , 



(60) 



for any four-dimensional tensor Z£. Using this notation, 
the Taylor expansion simplifies to: 

V 6 {n) °K^' = - V 1 QW (61) 

^{2n+l)\ v Z^( 2n )! vK ' 

n>0 v ; n>l y 1 



E 



L{2n) 



■o ( ™V„, 



which may be symbolically written as: 

— t= sinh \f£> °K$' = \ ( cosh \fd-\)8%. (62) 



To leading order in d, the derivative of the extrinsic cur- 
vature on the positive-tension brane is then: 



Kt'(y = 0) 
with 



■ tanh 




8» 



Mz>" :;i 



F^ = 65^ = 2d»dd v d - (ddf 5% (64) 
OFff = (ddf S» . (65) 

Therefore to leading order in d/L, the Weyl tensor on 
the positive-tension brane is: 



E i+)t ! = -~K<i'(y = 0) 



1 



D^D u d 
d ' 



D^D v d 



d 

1 

ML 
1 



(66) 



| ( M(tanhM_ tan MU 



, \ dd \ 

\dd\ v tanh — 

dL \L 



tan 



\dd\ 



with \dd\ = yJ-gW^dy.ddvd 
and Fj} = 2d> 1 dd v d - {ddf 5g, 

here the covariant derivative and the raising of indices 
are performed with respect to gffl = q^(y = 0). Note 
that, because of the modulus signs, we have 



(ddf = g i+) ^d^dd„d : 



(67) 



We may note that these functions are well defined only for 
\dd\ < ir. For \dd\ > ir, the inverse of the power expan- 
sion on the left hand side of I|61II62|I is ill defined. For the 
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purpose of this study we shall therefore restrict ourselves 
to the case where \dd\ < ir for which the expressions l)66|) 
is well defined. Whether or not the theory possesses ac- 
tual velocity-dependent divergences as \dd\ — > tt would 
be interesting to investigate with numerical simulations. 

All this analysis could be repeated using the negative- 
tension brane as reference brane instead. This would be 
equivalent to taking —d instead d all the way through. 
We would have then obtained the induced Weyl tensor 
on the negative-tension brane .E as 



E 



(-)*». 



D^D,.d 



d 

1 

+ 2dL 

+t— | tanh 
\dd\ \ 2 

1 J d 

H Ox - 

dL \L 



(68) 



\dd\ ( tanh 
\dd\ 



\dd\ 



where the covariant derivative and the raising of indices 
are now with respect to g^v — q^v{y =1)- This proce- 
dure has allowed us to derive the Weyl tensor on both 
the positive- and negative-tension branes in a covariant 
way. These expressions are exact in the close-brane limit. 
For simplicity we exposed in detail the case for which 
A 2 = g yy was assumed to be y-independent. It is unclear 
whether such a gauge can always be fixed, but our final 
answer is gauge-invariant and we show in Appendix I A 21 
that taking A to be y-dependent does not in fact affect 
the final result. For a discussion of this gauge issue see 

El 



IV. THEORY IN THE CLOSE BRANE LIMIT 

A. Modified Einstein equation 

Using the result of the previous Section, in the close 
brane limit, the positive-tension brane geometry will sat- 
isfy the following modified Einstein equation: 



w (IV 



_£(+) 

D^D v d 



1 

dL 



|Sd|tanhM 5 (+) 



(69) 



" w\ ( tanh 



\dd\ 

2 



tan 



\dd[\ 



dndd w d 



describing gravity coupled to a scalar field in a non-trivial 
way. This theory has some higher-derivative corrections 
as expected but, remarkably, they are both simple in 
form and entirely first order, involving only powers of 
first derivatives. This leads to the important result that 
the theory remains second order in derivatives even when 



correct to all orders in velocity. The initial data on a 
Cauchy surface that needs to be specified for the theory 
to be solved is the same as that needed in the low-energy 
effective theory. There will be no need to specify extra 
information or to consider the corrections to be small, 
a feature which will make the theory straightforward to 
solve. 

In order to solve this theory, we need to specify the 
equation of motion of the scalar field. So far the traceless 
property of the Weyl tensor has not been used, which 
gives rise to the modified Klein Gordon equation for the 
radion: 



Dd = 



\dd\ 
L 
1 



dd „ , dd 
tan - — - - 3 tanh - — - 
2 2 



(70) 



Using this equation of motion for the scalar field and the 
modified Einstein equation (|69|l . the geometry on the 
brane can be found in the close-brane limit. We may 
point out that only four-dimensional quantities are in- 
volved. Although this theory does not seem to be deriv- 
able from a four-dimensional action, it may be solved 
using standard four-dimensional methods. 



B. Low-energy limit 

In the low-energy limit, the four-dimensional effective 
theory presented in Section III B II is exact. For consis- 
tency we may check that, in the small velocity limit, we 
recover the small distance limit of this theory. To leading 
order in velocity, we may use tan ~ tanh ~ i \dd\ 
in H69fl and (JTftJ, giving rise to the Einstein equation and 
the equation of motion for the scalar field 



(+) _ D^Dyd 
" v d ' dL 



l -^dd v d- l -{ddfglt> 



Ud 



{ddf 
L 



which is precisely the leading order in d/L of the effective 
theory 1)251 126fl . The equation of motion for the scalar 
field is actually exact to all orders in d/L. 

We might think that the low-energy effective theory 
(|2*5)l could still give a correct answer provided that the 
relation between \P and d was modified. If that was the 
case, the relation between these two variables should in- 
clude some higher order velocity terms: W = ^(d,dd) 
such that, in the low- velocity limit, ^(d : dd) ~ e~ d l L . 
The theory should then include some terms of higher or- 
der in derivative of the form d a dD f2 D u d a d ^,dd- This is 
not compatible with the theory in (|69|l which only con- 
tains terms up to second order in derivatives. The close 
limit theory in Q69[) and the low-energy effective theory 
(|25|l are therefore genuinely different. 
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C. Conservation of Energy 



in velocities and in d/L, should be: 



In order for the theory to be consistent, the Bianchi 
identity on the brane needs to be preserved. This means 
that, to leading order in d/L, the traceless condition H70|) 
must imply the right hand side of (|69J) to be transverse. 
In other words, if we consider the right hand side of i|69|) 
to be the stress-energy tensor of the scalar field d, the 
Klein Gordon condition should impose that it is con- 
served, i.e. divergenceless. Using the equation of motion 
for the scalar field, we may rewrite l|69l) and (|70() in the 
form: 



E„ 



G (+) 
T 



DuDvd 
_d 

i r 



Ld 



f(z)giV+g(z)F^) (7i) 



D^Dyd + —gd^dd u d - dE^ u 

ad + ^ + 3 J--dE) 9 lV 



ad 



L 



(4f + 2z 2 g) +dE, 



(72) 



(73) 



\dd\, 



where for simplicity we used the notation z 
f = f{z A ) = § (tanh | — tan §) and g = g(z A ) = 
(tanh | + tan |) . / and g are viewed as functions of 
z 4 = (<9d) 4 since their Taylor expansions only contain 
powers of z 4 . Therefore, for example, 

d,f(z) = 2z 2 f(z)d,(z 2 ). (74) 

We have as well added the next-order correction dE^ LV 
where E^ v ~ d°. As we shall see, although E^ v does not 
contribute to the leading order of the Einstein equation, 
it does contribute to the leading order of the divergence 
of T^y. This is due to the fact that, although dE^ ~ d 
is negligible, its derivative is not (d^d) E£ ~ d°. We 
therefore need to consider its contribution as well. The 
transverse requirement will therefore impose a condition 
on the next to leading order in d. Using the result of 
Appendix [Bj the divergence of the stress-energy tensor 
as defined in lj?T|) is given by: 

dD^T}} = -^[gf + z 2 (f~z 2 g)(f + z 2 g')]d„d 



[EZ-E5t]d^d+—0{d/L) 



(75) 



= 0. 



For the stress-energy to be transverse to leading order in 
d we therefore have a constraint equation on the next-to- 
leading order contribution E^ v . Furthermore, from the 
low-energy effective theory (using l |25t ). to leading order 



L 



D^D v d + j- [ d^dd u d - ^ {ddf 4+) 



(76) 



~o([ddf,d/L 



A natural ansatz for Ea V is therefore: 



.(77) 



Using this ansatz in the constraint 175|) , we find the equa- 
tions for / and g: 

f~f + z 2 (g-g) - 4 (gf + z^fgi - z V) (78) 
+ 4, 2 (//'-zV). 

Since / and g are functions of z uniquely, their expres- 
sion is: 

/(z 4 ) = / + 4.g/ + 4z 4 (fg' - gf) (79) 
3(z 4 ) = g + iff - iz 4 gg'. (80) 

In the slow- velocity limit, we may check that / = 
0((ddf) and g = \ + 0({ddf) so we recover the re- 
sult from the low-energy effective theory ()76|l . 
Furthermore, using the expression H71|). the next to lead- 
ing order E^ (|77|) will be conformally related to the lead- 
ing order E^ v = j;E® only if / = / and g = g. But 
from the expressions (|79l I80|) , this can only be the case 
in the low- velocity limit where / = / = — 4p + 0(z 8 ) and 
g = g = 4 + 0(z A ). The higher order terms in velocities 
do not match and the conformal relation E^ = ^E® 
does not hold. 

This expression for E^ v with the relations (|79l I80|) is 
consistent with the close limit theory since it ensures that 
the stress-energy tensor is conserved to leading order in 
d/L. However we may emphasise that this is not the 
only possible answer and relies strongly on the ansatz 
(f77|l . Since we only want to focus in the leading order in 
d/L of the theory, the exact expression of E^ u will not be 
relevant, it is enough to show that it is possible to intro- 
duce its contribution in such a way that the transverse 
condition is preserved. If we wanted to go further in the 
study of the next to leading order term of the theory, 
we should check that this expression of E^ u is consistent 
with its exact expression that can be derived in the back- 
ground using the SAdS five-dimensional geometry. This 
is beyond the purpose of this work which aims to focus 
on the close brane limit and to neglect any terms which 
are not of leading order in the expansion. 



D. Consistency check for the background 

Having a theory which should be exact in the close 
brane limit, we may now check that it correctly repro- 
duces the exact background behaviour obtained in J7J 
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1201 • Taking the positive-tension brane to be endowed 
with a flat FRW induced metric with scale factor a + (t), 
the Einstein equations reduce to 



3H 



— I tan — 

d Ld 2 



(81) 



2H+ + 3Hl = H + - d - ± tanh | , (82) 



and the Klein Gordan equation l|7()J) gives 



d + 3H + d = — ( 3 tanh — — tan ^ 



Combining (|8*TJ) and gives 



Hi = -( I tanh- -flu 
+ d l L 2 4 



(83) 



(84) 



Also we have immediately that a+(t) is linear in confor- 
mal time (which is merely a consequence of the induced 
metric being FRW with empty branes). As with the anal- 
ogous low-energy result (I29|) . (|84|l has two solutions, one 
of which has d = at the collision and for which both 
d and H+ are finite and non-vanishing. Again, it is the 
latter we consider since this is the one corresponding to 
branes moving in opposite directions. Therefore for this 
study both H+ and a+ are finite in the limit d — ► 0, 
giving the leading order behaviour of H + from l|84|) as 



H+ = Z tanh 2 + V / 



(85) 



Replacing d with — d will give the corresponding result 
for the negative-tension brane: 



1 , d 
— tanh — 
L 2 



O 



d 



(86) 



The modified Friedmann equations on both branes (|85l 
I86|l are precisely the exact result (|20[1 we obtained from 
solving the five-dimensional equations in Section III Al 
This is a non-trivial check on the validity of this close- 
brane four-dimensional theory as it reproduces the right 
result to all order in velocities in the small distance limit 
and not only to leading order in d as was the case for the 
low-energy effective theory. 



V. EXTENSIONS 

Introduction of a potential to low-energy 
effective theory 



The low-energy presented in Section III B II was ex- 
pressed in the positive-tension brane frame. In order to 
add a potential to this theory, it is more natural to work 



in the Einstein frame. The Einstein frame is related to 
the brane frame by the conformal transformation: 



« = (cosh(>/V6))V, (87) 



<J 



g$ = (-sinhU/Ve) 



where any "bar" quantity designates a quantity with re- 
spect to the Einstein frame. The conformal scale fac- 
tor in H25fl is related to the scalar field cj> < by 
* = er d l L = - tanh (0/V6)- In that frame, the four- 
dimensional effective theory can be derived from the ef- 
fective action of a scalar field minimally coupled to grav- 
ity: 



(89) 



giving rise to the equation of motion for the scalar field: 
□ = 0, Cf. 0, E3| . 

It is possible to modify this effective theory by adding 
a potential by hand in 



8 = 



-9 



-R 



(90) 



In that case the equation of motion for the scalar field is 
modified to: 

U4> = u'{4>), (9i) 

and the Einstein equation in Einstein frame becomes: 

G>„ = d^d v( /> - Q(d0) 2 + Uj g^. (92) 

We can now go back to the positive-tension brane 
frame by performing the conformal transformation 18 7|) 

with (cosh (0/\/6)) 2 — 4i- To linear order in d/L, the 
theory on the brane is modified to: 



= -[D^D v d~UdgW 

+ j [d ll dd u d+- {ddY g^ v 



(93) 



J -< + >l-v(d) fl w 



□ d = 



L 



{dV'{d)-2V{d)) 



(94) 



where all covariant derivatives and index raising is taken 
with respect to gffl and for simplicity we used the po- 
tential V(d) defined as: 

V(d) =d(l 



* 2) U ~ — U. 
L 



(95) 



The modification of the equation of motion of the scalar 
field (|94|l ensures the right hand side of the Einstein equa- 
tion l|93|l to stay conserved even after the addition of a 
potential. The same procedure can now be applied to the 
theory in the close brane limit. 
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B. Addition of a potential in the close brane theory 

We consider the addition of a potential in the close 
brane theory by modifying the stress-energy tensor in 
the Einstein equation (|72|) : 



T — 



nd 



L 



g d^dd v d — dE^ 



(96) 



^^l-dE + Vid))^ 



For this stress-energy to remain conserved after the addi- 
tion of the potential, the equation of motion of the scalar 
field (|75)l has to be modified as well: 



□ d = -i(2z 2 . 9 + 4/) 



■dE + W, 



(97) 



where the correction term W is a functional of the poten- 
tial V(d) to be determined. The modified Klein Gordon 
equation l|97l) should be consistent with the conservation 
of the modified stress-energy tensor in (|96[1 . The diver- 
gence of T^j, is calculated in the appendix [B] using the 
constraint (|76|) for E^: 



dD^Tl 1 = d u d -V'(d) + 



3W + W 2g 



2d 



12 



+ T z 2 (/' + z 2 g') (V + W) 



+ ^-(3W + 2V) 
±0(d/L). (98) 



To lowest order in d/L we therefore need to modify the 
equation for the d with a term W: 



W = - (dV'(d) - 2V{d)) . 



(99) 



G 
I 



G 
I 



(100) 



The analysis is only slightly modified from the above and 
we shall only sketch it here for brevity. Using 143(1 . 152(1 
is modified to 

Y,h QK " (n) = \ {a++a - )5 ^ (101) 

n>l 

which gives rise to a more complicated version of ((63(1 : 



A;; (. v = in = + 



<T_|_ + O — 



O tanh 



-St 1 ,. 



6 sinhVO 
The analogous result to i(66|) is then 



(102) 



d vK!J ' d Z \ L 12 

D^D,,d 



d 

±(L + Z± 
2d \£ 6 



(103) 



j tanh M_ tan M),v: 



|9d| ( cosech |<9d| + cosec \dd\ ) <5^ 



(cosech \dd\ - cosec |dd| ) F? 



This is precisely the same term which had to be intro- 
duced in the equation (|9"4"}l for d in the low-energy the- 
ory. This procedure is therefore completely consistent 
with the low-energy theory and will give the same result 
to leading order in (dd). It is possible to modify the 
close brane theory by introducing a potential, in such a 
way that the conservation of energy and thus the Bianchi 
identity remain unaffected. 



C. Detuned tensions 

A full analysis of how the argument of Section llVl would 
run in the presence of matter on the branes is beyond 
the scope of this paper. We can, however, consider the 
simple example of brane cosmological constants where 
the stress-energy tensors are proportional to the induced 
metric to get an idea for how the coupling to matter 
might look. Specifically, we allow now for a detuning of 
the tensions to 



\dd\ 
3 U 12 



dL \L 



The Einstein and Klein-Gordon equations are then mod- 
ified to include complex couplings of the tensions to both 
the radion d(x) and powers of its first (but not higher) 
derivative. 



G (+) 



.#) - o(±) 



l 







a- 





\dd\ tanh 



2 " 



J_ hanh iyi + tan M ) n,,diU 



(<r + +a_) 



(104) 



\dd\ cosech \dd\ g^) 



+ 



1 

W\ 

2 
3 



(cosech \dd\ — cosec \dd\ ^ d^ddyd 



L 



12 



(105) 
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, dd \ dd \ 
3 tanh - — - - tan - — - 

2 2 



6 



\dd\ 3cosech \dd\ + cosec \dd\ 



_4d/^ o+\ 1 d 
In the low-energy limit, these reduce to 



(106) 



(a++a-) (+) , 1 
6d ^ 



□d 



(9dT 



cr+ +cr- 



(108) 



which agree with the standard results of the low-energy 
effective theory in the presence of detuned tensions |27j . 



VI. CONCLUSION 

In the first part of this work we pointed out the dis- 
crepancy between predictions from the effective four- 
dimensional low-energy theory and the exact five- 
dimensional results. The difference has been established 
in the regime where the branes were close and all through 
this paper we worked to leading order in d/L. 

In order to go beyond the low-energy effective theory, 
we established a formalism to find the electric part 
of the five-dimensional Weyl tensor on both branes. This 
tensor represents the only quantity which is unknown 
from a four-dimensional point of view as it encodes in- 
formation about the bulk geometry. Finding its expres- 
sion on the brane is therefore the key element in order to 
study the geometry of a brane within an extra dimension. 

Using this formalism in the small distance limit, we 
found an exact expression for the Weyl tensor on each 
brane, valid at leading order in d/L but at all orders in ve- 
locities, or for any energy scale. We were therefore able to 
modify the regime of validity of the effective theory from 
a low-energy regime valid independently of the branes 
distance to a regime valid at all energies for close branes. 
This regime of validity is relevant for cosmology since it 
represents one of the main focus of present braneworld 
models. Understanding the behaviour of branes just be- 
fore or just after a collision is indeed crucial if our Uni- 
verse emerged from a brane collision. If the Big Bang 
were initiated by such a collision, it is consistent to as- 
sume that the large scale structure was produced in a 
regime where both branes were close. Even if this regime 
is not valid after a while, effects that originated dur- 
ing this period are unlikely to be eliminated once the 
branes are far apart and the fifth dimension has opened 
up, leaving the possibility of exciting the Kaluza Klein 



modes. We therefore believe that understanding the con- 
sequences of cosmologies in this limit will allow us to 
study the viability of such scenario s by comparing them 
with observational data [TtL l2ll IM Ha IHol iHlf . 

In this paper we have established a theory that allows 
us to study such scenarios and have checked the consis- 
tency of its predictions for homogeneous and isotropic 
geometries. We argue that this theory will be remark- 
ably straight-forward to use since it includes only four- 
dimensional quantities and is effectively second order in 
derivatives, the only higher-order corrections coming as 
powers of first derivatives. This feature will facilitate any 
comparison with other four-dimensional theories. It is 
different from other higher-derivative theories pH 153. l53| 
in that it is not derived using the assumption that it can 
be derived from an action, and it is different from a the- 
ory relying on purely string effects EE EH • Already 
for the background, an interesting result can be pointed 
out: when the branes are empty the Hubble parameter 
on each brane is bounded by L 2 H 2 < 1, which could 
not have been derived from the low-energy effective the- 
ory. This does not hold anymore when some matter is 
introduced on the branes or if a potential is added for 
the radion. Another interesting result of the theory is 
that the expansion seems to break down when the veloc- 
ities are of order \dd\ ~ it. This is a feature which has 
never been pointed out before and we believe it suggests 
the presence of interesting physics which might be worth 
studying. 

To extend the theory, we considered the addition of 
a potential for the radion. This will be interesting if we 
want to study curvature perturbations as the origin of the 
large scale structure. In order to begin to understand the 
way matter on the brane would be coupled in this theory 
we have extended the model to the case where the brane 
tensions are detuned. Both these extensions will be useful 
for further progress in the understanding of braneworld 
cosmology after or before a collision. 
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APPENDIX A: LEADING ORDER DERIVATIVE 
OF THE EXTRINSIC CURVATURE 

1. Gauge where A is y-independent 

The purpose of this Section is to formulate the recur- 
sive expression for the derivatives of the extrinsic cur- 
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vature on the brane in the regime where the branes are 
close d<^ L. 

We recall from Section UlI C 31 that on the positive ten- 
sion brane, for y = 0, we have i^( m ) ~ d° for m = 0, 1, 2. 
Furthermore, from (|47jl . at y = 0, we have r' ~ dd ~ d° . 
We may as well calculate its second derivative: 



r M V'(y) = (d^');M + (d^');^( rf ^) ;a ( A1 ) 
-d<? a/ V 



1 I/p "p £± pLV IV p j 



r CT ft'' 3 

1 /3p /Y iy 



On the brane the last line is of order (d/L) 2 and the 
two middle ones are of order d/L. They are therefore 
negligible. The first line gives a contribution of order 
r"(0) ~ 8dK'(0) ~ d°: 

rr (Q) = ^dA-s'(o) + 0„d^'(O) - d a dK'^(o),m 

where terms of order 4 have been omitted. 

In what follows we will use a symbolic notation omit- 
ting any indices or coefficients. We write symbolically, 



lap 



9, K% 



K, d a = d and cL 



In particular, 



we have: q' = dqK. Using this notation and expressing 
E£ in terms of K'$ in l)48[l . we may symbolically write 
an equation for K" p : 



K" 



d (d 2 + 
+dKK' 



Td + dT + T 2 ) q' + dK 3 + dK (A3) 
+ q dd (T' + d TK + d dK) + d 2 d q', 



where again by d we designate a normal derivative (as 
opposed to covariant derivative) with respect to any co- 
ordinate x* 1 along the four-dimensional y = const hyper- 
surface. This expression is true for any y. The lead- 
ing order in d/L in this expression comes from the term 
qddV ~ d°. All the other terms are of order d. 
For n = 2, we therefore have on the brane: 



(0) - r (m ) (0) ~ d° V m < n. 



(A4) 



Now let us assume that this relation holds for a given 
7i = l + l. In particular, r( i+1 )(0) - d° and 



g( m + x )(G) ~ dK^d V < m < I + 1, 
q(0) ~ d°. 

For the next order in n + 1, we have: 
K {n+i) _ g{l)x" 

d (d 2 + Td + dT + T 2 ) q' 
+ dK 3 + dK + dKK' + d 2 dq' 
+ qdd (T' + dTK + ddK) 



(A5) 
(A6) 

(A7) 
(A8) 



with 



d^dd 2 q' 



dd 2 



1=0 



(0) 



(A9) 



Similarly we may check that 

flW [d (ra + ar + r 2 ) g' + d^ 3 + d/f] 

<9« [djjfif' + q dd (dTK + ddK) + d 2 dq'} 
and 



y=o 



d (A10) 
d, (All) 



0< z) [qT'dd] = dd V C; m Q (i - m) (0)r( /+1 )(0), (A12) 
w=0 m=0 

with ( ? ('- rn >(0)r( w + 1 >(0) - d 1 if m < I and 9 r(' +1 )(0) ~ 
<i so that: 



3^ [qT'dd] 



d U . 



v=o 



(A13) 



We have therefore shown that K^- n+1 \0) ~ d° and its 
leading contribution comes from the derivative of the 
Christoffel symbol uniquely. 

We want now to show that the same is true for the 
Christoffel symbol. We have K (m ^(0) ~ d° for any m < 
n + 1 and r( m )(0) ~ d° for any m < n. Using the relation 
(|47l) for T£j, we have: 

r;, (n+1) (y) = (d^ (n) ),p + idK a W\ v - q w {dK p v {n) y a 



T« p K + TZ p K^-2T(> v K<* p (A14) 



71-1 



m— 

where d^ n) (TK) ~ d. ^" _m) (rtpp) ~ d for n > m 
so the last sum goes as d. Finally the first term goes as: 
(dKZ (n) ),n = d tfl Kl [n) + 0(d). The Christoffel symbol 
therefore goes as r( n+1 ) ~ d° for all n and its leading 
contribution is: 

r;„ (n+1) (0) = d. M KZ {n) + d >u K a ^ - d' a K$ (A15) 

where terms of order 0(d 1 ) have been omitted. We can 
therefore finally conclude that 

#(») ^ r (n) - d° for any n > 0, (A16) 

and the leading contribution for the extrinsic curvature 
comes from the derivative of the Christoffel symbol only. 
Indeed for n — 2, we have: 



Kt"{Q)=q^Tfl d a d + 



L 



(A17) 



and similarly for any n > 2, the leading term in the 
extrinsic curvature comes uniquely from the derivative 
of the Christoffel symbol: 



K?S n) 



y=o 



q 0v T^ n - 1) d a d+O[y) (A18) 



= U^K /,^-^ +d^K atl{ - n - 2 ^ (A19) 

- d- a K^ n - 2 ^ d a d + o(^\ . 
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We therefore have a recursive expression for the leading 
order in d/L of the normal derivatives of the extrinsic 
curvature that we can use in order to formally sum the 
Taylor expansion (|45|) or 152fl . 



2. Gauge where A depends on y 

The previous result has been derived assuming that 
A in l|39[) was independent of y. It is possible to show 
that such a dependence does not affect the final result. 
In order to show this we will dissociate the derivatives 
of A from the other derivatives. Then by summing over 
the derivatives of A, we will recover the quantity d. We 
denote by S A = A' (y)8 A + A" (y)8 A , + d a A' (y)d daA + ■■■ 
the y-derivative which acts exclusively on A. 

In order to differentiate any other quantity, we can 
apply the operator d y — d y — 6 A which has no effect on 
A. In particular we write: 



B y Q(y) = d y Q(y) = Q'(y), 



(A20) 



for any quantity Q — K 1 ^, E^,,q^ v ,T^ v . In the previous 
Section, A was y-independent and so the operator 5 A 
had no effect, we simply had d y = d y . Now this does 
not hold anymore and we need to include the effect of 
S A in the Taylor expansion l|45l) baring in mind that S A 
and d y do not commute. Indeed, a new A appears each 
time the operator 8 y is applied on a quantity Q since 
B y Q = AL n Q, where L„ = 9j is the Lie derivative along 
the normal direction. So S A Q — whereas 5 A d y Q = 
OA'(y) where O is some four-dimensional operator. If 
we consider for instance the y-derivative of the extrinsic 
curvature in 144|) . we have: 

S A B y K^(y) = 5 A Kt'{y) 

= -A'(y)E2(y) - D»D v A'(y) (A21) 

-A'(y)K^y)K^y) + A'(y)^, 

In the Taylor expansion, (|45|l . there will a be a term, 
which we denote by JC^iO), that has only the first deriva- 
tive along d y and all the other ones along S A : 

= E ^4 ro_1) (A22) 



m>l 



We may write 

S^dyK^y) = 6^ m - x \ (A23) 

where he operator can be read off from (|A21Jl . Using 
this relation, we have: 

^-^ mi 

rn>\ 

= 6? [ V+ A(y')dy'. (A24) 



In particular the term ^^(0) that contributes to the 
Taylor expansion (|45(l is: 

£it(0) = OZ(y = 0) [ A(y)dy (A25) 
Jo 

= -dEZ-D^D u d\ y=0 , (A26) 

which is precisely the first term (|53|l that was contribut- 
ing in the Taylor expansion when A was assumed to be 
y-independent. 

It what follow we shall see that the Taylor expansion 
can be expressed as a sum of K, n which expression is 
precisely the same as K ^ when A is assumed to be y- 
independent. 

We shall work in a symbolic way, omitting any indices. 
First we notice that the Taylor expansion (|45l 15211 can be 
written in the form: 



1 



^ T7.I 

n>l 



-X n (0) = 0, 



with 



- y - 

mi,- jn n >{) 



(A27) 



(A28) 



where we wrote /„ = {mi + • • • + m 



n) and 

c(«h) • 



Pr'9v&A% ■ ■ ■ d v S A 'd y K. (A29) 

Recalling that each time K is differentiated with respect 
to y, a new power of A appears, so that 



AOr, 



(A30) 



where each operator Oi depends on y but includes only 
derivatives along the y = const hypersurface and quanti- 
ties Q = K,E, q, r. Using this notation, T n is: 

(A31) 



T, 



■mi , ■ ■ ■ ,m n 



Or, 



AO n -! A---6 2 (AOxA^ 



(m 2 ) 



(m„_i)' 



(m n ) 



where the operators Oi should be treated as independent 
variable of y so that for instance 

( J 40i^( mi ') ( " l2) = Y, ck m 2 A (k) d 1 A (mi+m2 - k K(A32) 

When all the sums are performed, we obtain: 
1 

n! 



—,fcn{y) 



(A33) 



O l 



A0 2 A 



AO n 



Ady )dy)dy 



So on the brane, we can express IC n in the form: 



dy. 



(A34) 



/C n (0) = n!f/i(l) 
U m (y) = 6 m / A(y)U m+1 (y)dy, m<n (A35) 



U n (y) = O n / A(y)dy = O n d(y). 
Jo 



(A36) 
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So far this result is exact. When the branes are close, 
this result will simplify remarkably. First we point out 
that < d(y) < d < L for any < y < 1. This follows 
from the fact that the bulk geometry is completely reg- 
ular, so between the branes, A(y,x) > 0, which implies 
< f£ A(y,x)dy < d{x) < L for any < y < 1. This 
holds similarly for any multiple integral. The same sim- 
plifications as in the previous Section will therefore be 
valid here. For instance if in the previous Section we had 
Old ~ Dd, then this will remain true in this case as well: 
6 1 ^A(y,x)dy^D^A(y,x). 

In the following we shall label with a "bar" any quan- 
tity which was derived in the previous Section assuming 
A was independent of y. In the previous Section, we had 
A = d and so U m was simply: 

-1— m 



U m {y) 



V 



[n + 1 — m)! 



0„ 



dO m +i (d.»<5nd)] , (A37) 



so we had: 



X (n) (0) = 6i d6 2 (d---d6id) 



(A38) 



using the fact that at d <C L, the action of the operators 
Oi was considerably simplified and at leading order in 
d/L, they were equivalent to the overall operator O^: 



with 6( 2n+1 ) 

Q(2n) 



OW<f 

1 



(2ra 
1 

2d 



D! 



(-25(0) 



if(0)L> (2n) , 



(A39) 

D 2^ D (2n) ( A4Q ) 

(A41) 



where Z) is a derivative along the y — const hypersurfacc. 
The leading contribution was indeed: 



£(2»+l)( ) 

K {2n \0) 



K'(0) {ddf 
K(0) (dd) 



2 n 



(A42) 
(A43) 



Now we may go back to the situation where A has some 
y-dependence. Because < A(y,x)dy < d(x) <C L 
for any < y < 1 and similarly for any multiple integral 
in U m , when the operators Oi are applied on these in- 
tegrals we can reproduce step by step exactly the same 
procedure as we followed in the previous Section in or- 
der to keep only the leading order in d/L. In particular, 
the repeated action of each Oi on each multiple integral 
can be substituted, in the small d limit (which implies 
a small Jq ■ ■ ■ Jq Ady limit) by the action of an overall 
operator O^ . If we do so, the leading contribution can 
be expressed in the same way as in (|A39() : 



Ady 



with the operator O^ as given in (|A40IA41|) . 
multiple integral is simply: 

a( f V ---A r Ady\ = i 

d n 



(A44) 



Now the 



-7 Ady 



(A45) 



So the rest follows exactly as in the previous case. For 
any n, the leading contribution to IC n is: 



/c n (o) = 6W<f, 



(A46) 



exactly as in l|A39(l . So the leading order in d/L of 
the Taylor expansion (|A27(I is independent of the y- 
dependence of A, and in particular we obtain the same 
result whether A = d or any other function of y such 
that Jq 1 Ady — d. The close brane limit theory derived 
in Section IIII CI is therefore valid independently of the 
precise expression of A, as this should be since our result 
is gauge invariant. 



APPENDIX B: DIVERGENCE OF THE 
STRESS-ENERGY TENSOR 

In order to work out in detail the divergence of the 
stress-energy tensor, we will consider directly the situa- 
tion with a potential in section IV Bl We therefore con- 
sider the stress-energy tensor as given in (|96|) with the 
equation of motion for d given in (|97() . We can work in 
the situation where no potential is present as in section 
(|IV C|l by setting V — W = in the following. 

Using ((HE!, the divergence of the stress-energy tensor 
may be written as: 

d Dfj/Ty = -V'(d)d v d-T^d fl d + R'id fl d (Bl) 
-i {6z 2 f + Az 4 g' + 2. 9 ) d v z 2 + ^g'F^z 2 

+~9D»F!t - d^d [Eg ~ E5at] - dD^ [Eg -ESS], 

where again / and g are taken as function of z 4 . We 
stress that Etf ~ d a so dD^E^ will be neglected in what 
follows. 

The first line simplifies, using the relation 

T£ = Rt+-T5>" v 

= RZ-^(3W + 4V)6ft. (B2) 

^From the Einstein equation and the equation for d (|97|l . 
we have: 



d^z 2 = -2D^D v dd v d 

= -2 (-\{fg$ + gF^) + (v + W)gj+A d-d 

= ^(f~z 2 g)d^d-2(V + W)d^d, (B3) 

where terms of order d/L 2 have been omitted. 
Using this relation and 



£> M F£ = 2Udd v d = -- (2z 2 g + 4/ - LW) d v d, (B4) 

Lj 
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the expression i|Bl|) simplifies to: 



[Eg - ESg] M d 
12 



(B5) 



-vy ) + ^ + 2 -l {3W+2V) 



+ ^z 2 (f + z 2 9 ')(V + W) 

±o ( d - 



d u d 



In the absence of a potential, the divergence of Tjf sim- 
plifies to the first two lines. Both 



12 



[Eg - ESg] d^d =-j^[gf 



(B6) 



L 2 

+z 2 (/-z 2 3 ) (f' + z 2 g')]d v d 
and, in the presence of a potential, 

v w . !E_t^ + | pw + 2V0 (Br, 

need to hold independently, in the close brane limit. 



APPENDIX C: EVOLUTION OF THE WEYL 
TENSOR 

In this Section we derive the evolution equation 148(1 for 
Eg. Essentially the same results are quoted in |10ll34Ll57| 
but we repeat them here for completeness. We assume 
the general form l|39l) for the metric, n = A~ 1 d/dy is 
the unit normal to hypersurfaces of constant y and, for 
a purely four-dimensional tensor, the Lie derivative with 
respect to n is given by 

L„-n = A~ x d/dy 

The projection tensor h a t — g a b — n a rib then satisfies 
/itj = /ig = = 0, hg = Sg and so can be left implicit 
in the index convention. The four-dimensional and five- 
dimensional Riemann tensors are related by 

Rnvafi = ^ R/Mvap + ^K^ a K^ v (CI) 

The electric and magnetic parts of the bulk Riemann 
tensor are defined to be 



E ab = (5) R acbd n"n 



B abc = h a h f b n d ^R efcd (C2) 

and are purely four-dimensional. 144|l then follows iden- 
tically from a direct evaluation of the Riemann tensor 
for the metric l|39l) . Denoting by V a the five-dimensional 
covariant derivative, we find that 

1 ft \ \ 2 

+2K^B v]Pa -K^B^ p {CZ) 



J {5) R„ 5ac n c ) = Dji„ u -K?.B. 



A 



vf3a 



(C4) 



As a consequence of the five-dimensional Bianchi identi- 
ties, we find 

V [5 ( (5) iV ]Q X) = (V [5 n c ) (5) i?^ ]QC 

which, from 1C3(I and l(C4|l . gives 



tin-B^Q, — K^B^a — —d^AE^a — 2D[ fi E^ a 

^ (5) iW (C5) 



A 



Similarly V^'R MV ] a p = yields 

2 2 

= n ■ d (5) R^ af j - —B^ [a dp]A- —Bap^d^A 

+2^ 5 'R fll/p [ a Kt + 2D[^B\ a p\^ + 2K ll [pE a ] I/ 

— 2Kv[pEa]n 

(tcTTl and JSJl can then be used to write this purely in 
terms of four-dimensional quantities, giving 

2 - 

Ln-R Mra /3 = -2R^ p [ a K^ + —B^adflA (C6) 

2 - 

— 2£)[ /x B| a( a|„] + -^Ba^d^A 

2 2 
-■jD^diaAKflv + -^D v d [a AK 0]fl 

All the above is true identically off-shell, i.e. not using 
the bulk Einstein equations ^ 5%> R a b — —jji9ab These relate 
the projections of the Riemann tensor to those of the 
Weyl tensor, 



E p u — ^ Cp,cvdji n> — E^v 



L 2 



B 



fiver 



( 5 ) C n c = B 



and we can use l(41|l to write 



^-'n-BjjiV — ^nffiu R[iv ^ lif 

f = KK — K K p 

The first term on the RHS can be found with repeated 
applications of (|44|) whilst the second can be found from 
(|C6|) : the result is 

= D a B a i ytlu ) + {Kq^ v - Kfj, u ) 
+K af3 R pav(i + 3K?„E v)a - KE^ 



This can be written in terms of quantities which vanish 
on the brane where K pv cx q pv , namely the traceless part 
Kg = Kg — jKSg of the extrinsic curvature and the Weyl 
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tensor (which will vanish when contracted with K al3 ): which reduces to i|48|) using -B M „ P = 2D^K v ] iP . 
L n E^ = 2K^E lux -^KE IM ,-^K^E^q iaf (C7) 

+c IMUf3 K a P + 2k«k af3 k« - ^k a/3 k a ^k^ 
-^q^k aP k^k a p + D a B a(pu) 
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